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COMPENSATION OP LARGE CHANGES IN REACTIVITY 
BY MEANS OF CORE LATTICE DEFORMATION 

G.A,Bat, G.B. Muhina, D.M.Parfanovich, D. C,Klochkov, 
V.C.Muhortov, H.M.Truhachev, D. I.Sheffer. 

1 . 

resuSts' t ' tem P' fc t0 prolong core of power and research reac- 

tors in continuous increase in excess reactivity of clean cores. 
In light water-cooled and-moderated reactors which are more 
widely used this excess reactivity is about 0,15 and tends to 
increase with every new core developed. At the same time methods 
of excess reactivity compensation are developed and improved: 
burnable poisons are inserted, mechanical systems (control rods, 
fuel assemblies, and comg|nsating ^ * eras- rods moved in groups) 
become greater in number more complex. The present methods, how - 
ever, have some disadvantages, therefore they cannot be conside- 
red quite satisfactory. Thus, the use of the most peri'ect system 
of burnable poisons cannot ensure the operating excess reactivity 
required for compensating temperature effect, steady-state poiso- 
ning^ power changes and "boron disbalance”, and mechanical systems i 
for example, compensating dystems, give rise to large power dis- 
tortions or as systems of a great number of low capture rods with 
separate drives are too complex and bulky. In addition a system 
of rods and their drives result in increasing the reactor height 
required. 

are three possible ways of providing operating excess 

rea.Ci.:.v. :j: 

a) Usual absorbing control system mechanisms: rods, assemb— 
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lies, compensating systems and their various modifications. 

b) Changes in moderator properties during reactor operation, 
for example, the use of light and heavy water mixturos of various 
component concentrations, solution of absorbing materials in the 
moderator and their gradual extraction, 

c) Changes in core lattice spacing (or some other deforma- 
tion of the core). 

The use of conventional control mechanisms is easily combi- 
ned with changes in the moderator properties but it is rather 
difficult to achieve in the case of core deformations. Any meth- 
ods of providing operating excess reactivity permit the use of 
burnable poisons. 

The possibilities of traditional compensation systems are 
well studied and practically they always prove to be rather 
"rigid", i.e, they contain a small number of highly efficient 
control mechanisms. Changing the moderator properties during the 
reactor operation is an example of "soft" control giving rise to 
minimum power distortions. Control systems of this type begin to 
find use in some reactors* 

The specific features of a possible control system hy means 
of core lattice deformation are being discussed below. 

Let us consider a reactor core formed by a regular lattice 
of cylindrical fuel elements placed in a liquid moderators, say 
in H 2 0. Let the top ends of all element e fastened in some plane 
thus that any of them can turn independently around the fastening 
point. Let, further, these elements pass freely through the spacer 
plate, the Mies in it forming a grid similar to that of fastening 
points. The lattice and the form of the cylindrical reactor do 
not depend on the distance between the plane of joints (fastening 
points) and that of the spacer plate. Let the plate be rotated in 
the horisontal pi, ie around the reactor central axis. Then the 
core will have the form of a one-sheet hyperboloid of revolution 
and the fuel element lattice will be deformed thus that the sur- 
faces of the constant spacing will turn out to be, with a greater 
accuracy, ellipsoids of revolution having one focus with the 
hyperboloid that forms the side boundary. 

With i* lie given fuel elements and the angle of twisting the 
reactor form uependa on the distance between the spacer plate 
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and the plane of joints, therefore control may be also provided 
by changing thib distance but not with rotating alone. In prin- 
ciple, a case may be considered when the spacing of holes in the 
spacer plate is different from that of joints and the initial 
core form is not a cylinder but a truncated cone. However, for 
design considerations as well as for improvement of power dist- 
ribution the variants symmetrical or nearly symmetrical relative 
to the core central plane are of greater interest. 

The proposed schemes are not suitable for reactors formed 
of separate fuel elements cooled by water because of changes in 
coolant passage area during the core deformation. However, this 
is not the case in systems with fuel elements enclosed in chan- 
nels as interchannel water is only a moderator. 

The reactivity dependence on the core spacing in a cylind- 
rical reactor is not monotone, it has optimum. When the spacing 
increases reactivity rises first, and then falls. Therefore to 
control a reactor by means of the core deformation it is possible 
to use either the "left" branch of the curve corresponding to the 
core spacing less than optimum or the "right" one corresponding to 
the spacing greater than optimum. 

It is evident in advance that the maximum efficiency of cont- 
rol for the right branch is limited by nuclear safety considera- 
tions. Really, in a heterogeneous reactor increase in lattice 
spacing causes increase in positive component of temperature 
effect due to flattening the thermal neutron distribution over 
the cell. In a system with a great spacing lattice the tempera- 
ture effect as well as temperature coefficient of reactivity 
becomes positive near the operating point. Keeping the conside- 
ration pointed out in mind it is desirable to make initial lat- 
tice spacing not very great. 

The minimum spacing for the left branch is determined accor- 
ding to the channel diameter and the possible deformation effici- 
ency depends upon the difference between reactivity with optimum 
spacing and that with the channels touched. In the case of chan- 
nels containing little interchannel moderator or in an extreme 
cas-. -f separate fuel elements the power distribution along the 
recO ... axis may turn out to be unstable while twisting. While 
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decreasing the coolant passage area the small asymmetry of the 
lattice and the external form of the core relative to its central 
plane give rioe to a non-uniform neutron flux distribution. 

In a reactor with a deformated lattice it is rather diffi- 
cult to locate conventional quick-response mechanisms of scram 
and automatic control circuits. A possible design variant may be 
performed by removal of channels from the region around the reac- 
tor axis or by using a group of central channels "twisted" or 
brounght closer in advance. In addition, one should care not 
only for efficiency of quick-response control mechanisms but for 
prevention of great radial neutron distribution distortions. But 
the greatest difficulty encounted may be the necessity to provide 
for reliability of a great number of kinematic pairs and flexible 
elements ensuring the movement of channels. In essence, feasabi- 
lity of the compensation system under discussion depends on the 
possibility to provide for a safe operation of these pairs and 
elements in the reactor^wftach is not accessible while in opera- 
tion. 


To describe mathematically the control system by means of 
"twisting" the core lattice, let us consider a cylindrical reac- 
tor, where 2H is an effective height oi v he reactor and R is an 
effective radius. If Z = + Zo plane is turned to + (j anglts 
relative to the vertical axis the cylinder will change into r- k 
body limited by a hyperboloid of revolution and a surface that, 
with reasonable accuracy, may be taken for an ellipsoid. In this 
body cross-section by the plane passing through the axis Z, AD 
and BC sections have the form of ellipsoid arcs 

J± „ At _ j 

A ^ f-Z - 1 


and AB and DC sedtions are the arcs of the hyperbola having one 
focus with the ellipsoid 


* = / 
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z # x 


Let the focus distance for the both of the cirves be 2c. It 
is easy to state for simple geometric considerations that 

d, -a = £c*fco 
4 - & 


l f * fa / £ 


p . 

% " 

do = 

$9 = ff 


__ l/-/^ +J- 9 


( 



As 


4 J.n . J z^.6, 


Thus the parameters a, b, c of the twisted reactor are comp- 
letely defined provided the effective height and radius of the 
cylinder, the rotation angle (a) , and the distance between the 
twisting surfaces 2Zo, are given. The maximum angle of rotation 
is defined by the contact condition of the fuel elements with the 
diameter d: 


fasui CO 


D is the cell spacing in the cylindrical reactor [u) = 0). 

The boundary problem equation of finding the critical para- 
meter in one-group approximation will be witten as follows: 





as 


Q3 ^ 
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( 1 ) 
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y j / 'z - o 

<» (2) 

where r is the thermal neutron flux, £ effective reactor 

surface, 2 material buckling depending on the apace co-ordina- 
tes, neutron flux and integral flux. 

The nonlinear equation (1) may be considered quasi-statio- 

nary. 

Either the angle U) or reactivity jO may be considered 
the critical parameter of the fuel burn-up problem. Ae compared 
with the results of critical experiments, it is convenient to 
consider the effective radius R the critical parameter and to use 
in the calculation 36- vqlue versus the cell area S determined 
experimentally instead of x 2 value calculated on the basis of 
geometrical and nuclear lattice parameters. This eliminates pos- 
sible errors in calculating multiplication factors and this makes 
it possible to evaluate easily the errors due to approximation in 
describing a real reactor system by means of one-group model for 
a reactor of effective sizes. 

It is convenient to solve the boundary problem in degenerate 
ellipsoidal co-ordinates for an oblate ellipsoid of revolution: 

X — C (Ufd (j? 

y = c aA«y3 (3) 

^ = c c*r , 3 


The co-ordinate surfaces are oblate ellipsoids of rotation, 
- = const, one sheet hyberboloids of rotation, J = const, the 
planes tf = const. The Laplacian operator in this system of the 
co - ordinates has the form: 







f * d 

Y <^C ci j l YYYi 
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Let us confine ourselves to the case when the solution 
does not depend on the angle q> . 

Then the equation (1) will take the form: o ^ » 

iU i ' if r ' 


-+- cJ 


/; £ / />■ c 2- / £i } . /' V' 

e-/C j- 0y ^ r 


It is easy to become eure that 


/ y a 


£ =- ^ e Z o£<^C c 06 Crtjd, -t 

__ 6S-J J%, 

a 

where sin - —3“ 

To use finite differencial scheme it is convenient to int- 
roduce variables 

^ ** 1 - cos J3 
O * ah cK 

and to transform the Equation (1) and the condition (2) to the 
form 

f 7 iJ if ]< ( 1 h , 

.s/fU-v*]**-?.*. 0 

<P= O when ^ = 1-cosy^,, when y = sh^« -5- 

. - ■— s 0 when iT = 0 when /? * 0 (2 ) 

Qn ^ ( 

(IA case of numerical solution it is possible to confine oursel 
ves to the quater of the symmetrical figure)# 

The expression for S in new variables is 

a c 7 vi ^7 


(2 11 ) 


^ & 






' £ + CtjJs, 


The nonlinear equation (1”) with boundary conditions (2") 


_ n _ 
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was solved numerically on an electronic computer. The accuracy of 
the solution was compared with analytic results for special cases 
of as a function of the co-ordinates. 

3. 

The suitability of the mathematical model (1) and (2) was 
analysed by comparison of predicted data with the experimental 
results obtained for a set of fuel elements containing enriched 
to 7% U0g encased in aluminium cans. The outside diameter of the 
fuel rods is 11.4 mm, they are 1200 mm long. The view of two assem- 
blies (one of them "in cross-section") is shown in Pigs 1 and 2. 

The critical numbers of fuel elements N of those assemblies 
with different lattioe spacing and angles of "twisting": were com- 
pared: « 

v jf f&W/tO ~ P J 

i3 % 


the d value describing the reflector saving was considered the 
function only of the initial spacing D, 

The relation a^(S) was obtained as a result of experiments 
on cylindrical systems. The agreement between the predicted and 
experimental results was quite satisfy tv. The calculated neut- 
ron flux distributions along the axis of the assemblies (Pigs. 3, 4) 
are rather accurate too. They were compared with the distributions 
obtained by means of a semiconductor crystal surface-barrier coun- 
ter recording fission fragments from a U-235 target. The discre- 
pancy in the distributions was only for the closely-spaced assemb- 
lies twisted thus that the fuel elements nearly touched each other 
(Pig. 5). Obviously this discrepancy was coused by non-ideal geome- 
try of the experiments: the lower boundary of the core was flat 
and the upper one was too convex. A plexiglass cap used as the 
bottom reflector considerably reduced the neutron flux distribu- 
tion asymmetry (see curve 3 in Pig. 5) though its curvature was 
noticeably less than it was necessary for the core to have a 
regular form. 

To calculate the reactivity value compensated by "twisting" 

- 8 - 
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with a fixed number of fuel elemental in equation (1) was 
aubatituted for 

3^ 2 -9^ (1 -f ) “//M 2 

The critical parameter, reactivity p , versus 00 at two 
values of the initial spacing D is shown in Pig. 6. 

In the whole range of 0^ 0.09 the axial non-unifor- 

mity of power distribution differs from the initial value for not 
more than 20%, the deformation of a closely-spaced lattice impro- 
ving and that of a great spacing lattice worsening neutron flux 
distribution. The corresponding parameters of volumetric non-uni- 
formity are -v 30% and ^10%. It should be noted here too that 
the fuel burn-up process promotes the power distribution flatte- 
ning. 

A second set of experiments was made on channels in the form 
of a bundle of concentric tubular fuel elements in steel cans 
inserted one into the other. One of the assemblies is shown in 
Pig. 8. As there s no water in the channel, vs.D had no going 
down branch at small spacing and it was possible to use only the 
range D of 70 to -\/ 50 mm. The JO value as a function of CO 
was studied experimentally. An additional channel was inserted 
into the critical assembly and reactor period was measured. Accor- 
ding to generally used formula the period was transformed into 
A j? t a, change in reactivity in case of an additional charge. 

The A JO value as a function of N was integrated numerically 
with the upper limits of 170, 160... 90, and 85 channels; integ- 
ral values are shown in Pig. 7. Additional experiments were made 
to evaluate an operation of possible shut-down systems. In par- 
ticular, the relation between the angle of , ’twipting' , and the 
quantity of boric acid injected provided N = const = 152 was 
studied. H^BO^ concentration was 0.192 and 0.396 g/litre at 
OJ * 15° and 20° • These values are very far yet from the solubi- 
lity limit and this makes it possible to use a small volume of 
highly concentrated solution injected into the core as a safety 
means in case of emergency. 

The central channel waB removed from the same critical as- 
sembly and a steel sheath tube, 44 o.d, 36 i.d, was placed instead 

of it and then a bundle of boron rods, 8.5 mm in dia.each, was 

376 
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inserted into the sheath tube* The experimental results are lis- 
ted in Table I. 


It is obvious that such a bundle or a single rod of the cor- 
responding diameter may be used as an automatic regulator. Simul- 
taneously a programme for design of systems that are two-region 
in the radial direction was prepared* 

It is shown that the radial flux distribution in an assembly 
with a cavity formed instead of seven central channels removed is 
quite acceptable. Neutron flux peaking in the cavity influences 
favourably the efficiency of absorbing rods that are inserted 
into it, in addition the rods may be fabricated rather thick. 
Unmovable channels for absorbers may be twisted in advance. The 
water gap between the absorber zone and the core may give rise to 
increase in power of the inner ring formed by 12 channels though 
decrease in nuclear fuel concentration prevents the emergency con- 
ditions from ocurring. 

The reactor scheme under consideration becomes more complete 
if fuel channels are in an annular position instead of forming a 
hexahedral lattice. Each ring should contain the same number of 
channels as the corresponding row of the hexahedral lattice does, 
that is 6, 12, 18, etc. 


Ap value in terms of 10”- 3 


Table I. 


y - 2 CO 

i 20° 

! 

O 

o 



l 

! 

j 

45° 

Central channel removal 

8.6 

10.7 


11.7 

Steel sheath tube insertion 

0 

0 


0 

Insertion of boron rods bundle 

3.3 

4.8 


6.3 


The distance between the rings is chosen thus that every 
channel will have an equal water volume. For example, the exter- 
nal ring radius in the direction of the channel axis is 

L — < N “ T~U where N is a total number of cells, n a num- 
ber of cells in the external ring. It was show l experimentally 
that this deviation from the regular form of cells does not affect 
the reactivity. 

376 - 10 - 
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Pig. 3. Changes in neutron flux distribution along assembly 
axis with 3 cm initial spacing while "twisting". 

Pig. 4. Changes in neutron flux distribution along assembly 

axis with 1.7 am initial spacing while "twisting". 
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Pig. 5. Neutron distribution instability at low and great 
D values 

1. Calculated curve for symmetrical assembly. 

2, Plat lower boundary of core. 

3.. Partially symmetrized assembly. 
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